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We present a model, motivated by the criterion of reality put forward by Einstein, Podolsky, 
and Rosen and supplemented by classical communication, which correctly reproduces the quantum- 
mechanical predictions for measurements of all products of Pauli operators on an n-qubit GHZ state 
(or "cat state"). The n — 2 bits employed by our model are shown to be optimal for the allowed 
set of measurements, demonstrating that the required communication overhead scales linearly with 
n. We formulate a connection between the generation of the local vàlues utilized by our model and 
the stabilizer formalism, which leads us to conjecture that a generalization of this method will shed 
light on the content of the Gottesman-Knill theorem. 

PACS numbers: 03.65.Ud, 03.67.Lx, 03.67.-a 
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I. INTRODUCTION 

Bell's theorem p| codifies the observation that entan- 
gled quantum-mechanical systems exhibit stronger cor- 
relations than are achievable within any local hidden- 
variable (LHV) model. Beyond philosophical implica- 
tions, the ability to operate outside the constraints im- 
posed by local realism serves as a resource for informa- 
tion processing tasks such as communication |2| , compu- 
tation [3J, and cryptography |^|. 

The violation of Bell-type inequalities demonstrates 
the in-principle failure of LHV models to account for all 
of the predictions of quantum mechanics. One approach 
to quantifying the observed difference between classically 
correlated systems and entangled states is to translate a 
quantum protocol involving entanglement into an equiv- 
alent protocol that utilizes only classical resources, e.g., 
the shared randomness of LHV's and ordinary classical 
communication. Toner and Bacon p| showed that the 
quantum correlations arising from local projective mea- 
surements on a maximally entangled state of two qubits 
can be simulated exactly using a LHV model augmentcd 
by just a single bit of classical communication. Piro- 
nio took this analysis a step further, showing that the 
amount of violation of a Bell inequality imposes a lower 
bound on the average communication needed to repro- 
duce the quantum-mechanical correlations. 

The original Bell-type inequalities Q, Q were for- 
mulated for pairs of qubits. Greenberger, Horne, and 
Zeilinger introduced a qualitatively stronger test of lo- 
cal realism, based on a three-qubit state, ^3) = (|000) + 
| lll>)/V / 2, which is now called the GHZ (or "cat" ) state. 
Here |0)(|1)) represents the eigenvector of the Pauli Z op- 
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erator with eigenvalue +1(— 1). GHZ correlations have 
been expcrimentally dcmonstrated in entangled three- 
photon systems and shown to be useful for performing 
information-theoretic tasks such as entanglement broad- 
casting [Hland quantum secret sharing 11]. 

Mermin |12| introduced a simple argument that shows 
how correlations between Pauli operators measured on a 
GHZ state violate local realism. We briefly review Mer- 
min's argument in Sec.[nj Mermin's formulation is based 
on the Einstein-Podolsky-Rosen (EPR) 0] reality crite- 
rion: "If, without in any way disturbing a system, we 
can predict with certainty . . . the value of a physical 
quantity, then there exists an element of physical reality 
corresponding to this physical quantity." This criterion 
is meant to capture what it means for a physical system 
to 'possess' a certain property. 

Taking the EPR concept of an element of reality as our 
starting point, we formulate a LHV model for the n-qubit 
GHZ state. By itself, the model is inadequate. It cannot 
give the correct quantum-mechanical predictions for mea- 
surements of arbitrary products of Pauli operators and 
correlations among such measurements, as is clear from 
Mermin's argument and its generalization to n qubits. 
Nonetheless, as we show in Sec. IIIII when the model is 
augmented by n — 2 bits of classical communication, it 
does reproduce all the quantum-mechanical predictions 
for measurements of Pauli products and their correla- 
tions. We go on to prové in Sec. II VI that this amount of 
classical communication is optimal for the allowed set of 
measurements, i.e., for measurements of Pauli products. 

In Sec. we demonstrate that our model arises nat- 
urally from a LHV simulation of a quantum circuit that 
creates the n-qubit GHZ state. The quantum circuit 
consists of an initial Hadamard gate H followed by a 
sequence of controlled-NOT (C-NOT) gates. It is a spe- 
cial case of a general class of quantum circuits identificd 
by the Gottesman-Knill (GK) theorem 0. The GK cir- 
cuits are those composed of qubits (i) initially prepared 
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in the state |00...0), (ii) acted upon by gates in the 
Clifford group, which isgenerated by H, 90° rotation 
about Z, and C-NOT [ÏJ], and (iii) subjected to mea- 
surements of products of Pauli operators. These circuits 
are capable of generating globally entangled states, such 
as the GHZ state, but their evolution can nevertheless be 
simulated in Q( n 2 / logn) operations on a classical com- 
puter [3L ÍT^. rÏ5|. We return to the question of GK simu- 
lations in Sec. IVII comparing and contrasting them with 
our simulation of the creation of a GHZ state and spec- 
ulating on how our results might impact understanding 
of the GK theorem. 



II. GHZ CORRELATIONS 

Mermin's three-qubit GHZ argument can be summa- 
rized as follows. The three-qubit GHZ state ^3) is 
uniquely specified as the simultaneous +1 eigenstate of a 
complete set of commuting Pauli products, one choice for 
which is (-XYY,-YXY,-YYX), where the ordering 
in the product specifies which qubit the Pauli operator 
applies to. In the language of the stabilizer formalism 
0, uM i the three commuting operators are referred to as 
stabilizer generators of \4>s). The stabilizer generators 
give the definite outcome +1 when measured, implying 
that a measurement of two of the Pauli operators in a 
generator can be used to predict the result of a mea- 
surement of the third with certainty. Thus, according 
to the EPR reality criterion, we should associate a local 
element of reality, having value +1 or —1, with the X 
and Y Pauli operators of each qubit. Letting Xj and yj 
denote the vàlues of these six elements of reality, where 
j labels the qubit, the stabilizer generators require that 
Xiy2V3 = U1X2V3 = V1V2XZ = -1. Multiplying these 
three quantities together gives X1X2X3 = —1, showing 
that the model predicts the result —1 with certainty for 
a measurement of XXX. Because of the anticommuta- 
tivity of the Pauli operators, however, the product of the 
stabilizer generators is +XXX, showing that quantum 
mechanics predicts the result +1 for this measurement 
with certainty. Mermin's GHZ argument demonstrates 
the incompatibility of quantum theory with local real- 
ism. 

The n-qubit GHZ state, \ip n ) = (|00...0) + 
|11 . . . l))/\/2, is specified by n stabilizer generators 
(X®™, ZZI® (n ~ 2 \ZIZI® (n - 3 \ . . . , Z/®("- 2 )Z), where 
/ is the identity operator. The full stabilizer group 0, 
generated by these generators, consists of the 2™ com- 
muting Pauli products of which \ip n ) is a +1 eigenstate; 
it contains Pauli products that have (i) only I's and an 
even number of Z's and (ii) only X's and an even num- 
ber of Y's, with an overall minus sign if the number of 
Y's is not a múltiple of 4. Of the 2x4™ Pauli products 
(including a ± in front of the product), 2™ are members 
of the stabilizer group, 2™ are negatives of the stabilizer- 
group elements and thus yield —1 with certainty when 
measured, and all the rest return ±1 with equal proba- 



X Y Z 

qubit 1 R 2 R3---Rn iRiR 2 ---Rn Ri 

qubit 2 R 2 ÍR1R2 Ri 

qubit 3 i?3 ÍR1R3 Ri 

qubit n R n iRiR n Ri 



TABLE I: Table of LHV's associated with an n-qubit GHZ 
state. Each row corresponds to a qubit, and each column 
to a measurement. The quantities Rj denote classical ran- 
dom variables that return ±1 with equal probability. The 
origin and meaning of the subscripts j becomes clear when 
we consider the creation of a GHZ state in Sec.[V] The out- 
come predicted for a joint measurement of a Pauli product 
is obtained by multiplying the corresponding table entries for 
each qubit (using 1 for unmeasured qubits, i.e., for an iden- 
tity operator appearing in the Pauli product) and discarding 
any factor of i in the final product. For example, for a joint 
measurement of XYY on the (n — 3)-qubit GHZ state, our 
model predicts the outcome (R2R3)(iRiR2)(iRiR3) = — 1, 
in agreement with quantum mechanics. Here we have used 
the fact that iï| = 1. Similarly, for a measurement of IYZ, 
the product of table entries is (iRiR2)(Ri) = 1R2; with the 
i discarded, the predicted outcome is the random result R2, 
again in accord with quantum mechanics. The use of i does 
not mean that the results of Y measurements are imaginary; 
rather the i is a "flag" that telis us how to combine Y vàl- 
ues for diflerent qubits in a joint measurement. Although 
the LHV table might seem not to respect the qubit-exchange 
symmetry of the GHZ state, one easily sees that it does by 
defining R' 2 — R2 ■ • ■ Rn, which exchanges the roles of the first 
and second qubits in the table. 



bility. 

Mermin's argument generalizes straightforwardly to 
\ipn) (our proof of optimality in Sec. IIVI can be viewed 
as just such a generalization) and shows that no lo- 
cal realistic model can correctly predict the outeomes 
of all measurements of products of Pauli operators per- 
formed on \tp n ) and correlations among such measure- 
ments. We now present a classical-communication- 
assisted LHV model that does yield all of the correct 
quantum-mechanical predictions. 



III. COMMUNICATION-ASSISTED LOCAL 
MODEL OF GHZ CORRELATIONS 

Our LHV model is specified in Table ^ which lists lo- 
cal realistic vàlues for the X, Y, and Z Pauli operators 
of each qubit. The caption describes how to determine 
the predicted outcome for a measurement of any Pauli 
product by multiplying the appropriate table entries and 
discarding any factor of i in the final product. The use of 
the imaginary phase i in the Y column, apparently just 
a curiosity, actually plays a crucial role. It reconcilcs 
some of the conflicting predictions of commuting LHV's 
and anticommuting Pauli operators, which form the basis 
of Mermin's GHZ argument. More precisely, the multi- 
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plicative àlgebra of these phases provides a concise rep- 
resentation of the n — 2 bits of classical communication 
required to ensure that our LHV model yields all of the 
correct quantum-mechanical predictions. 

To show that Table [ï]gi ves correct predictions for mea- 
surements of Pauli products, we consider those measure- 
ments for which the table predicts a definite outcome. 
Suppose first that a Pauli product contains no X's or 
y's, but consists solely of I's and Z's. Then it is clear 
that the table predicts certainty, with the outcome being 
+1, if and only if the number of Z's in the product is even. 
Suppose now that the product has an X or a Y in the 
first position. Then it is apparent that to avoid a random 
variable in the overall product, all of the other elements 
in the product must be X's or y's and the number of y's 
must be even; the outcome is +1 if the number of y's is 
a múltiple of 4 and —1 otherwise. Finally, suppose the 
Pauli product has an X or a Y in a position other than 
the first. Then the only way to avoid a random variable 
in the overall product is to have an X or a Y in the first 
position, and we proceed as before. These considerations 
show that our model predicts a definite outcome, with 
the correct sign, for precisely those Pauli products that 
are in the stabilizer group (or their negatives), includ- 
ing as a special case the observables forming the basis of 
Mermin's GHZ argument. Likewise, the model correctly 
predicts a random result for all other Pauli products. 

Our model correctly predicts the outcomes for all mea- 
surements of Pauli products, including single-qubit mea- 
surements. It fails, however, in some of its predictions for 
correlations between single-qubit measurements. To be 
correct, the model would have to reproduce all these cor- 
relations for all sets of single-qubit measurements. The 
model fails because products of single-qubit measurement 
results predicted by the model are not always equal to the 
corresponding joint measurement results. This inconsis- 
tency is a direct consequence of the rule that discards i 
from a calculated measurement outcome. As an exam- 
ple, consider the single-qubit measurements XII, IYI, 
and IIY on a three-qubit GHZ state. The product of the 
single-qubit measurement results, R2R3 for XII, R1R2 
for IYI, and R1R3 for IIY is +1, which is inconsistent 
with the prediction of the model and of quantum me- 
chanics for a joint measurement of the observable XYY. 
From the perspective of the model, classical communi- 
cation between qubits, an obviously nonlocal element, 
is necessary precisely to ensure the consistency of joint 
measurement predictions with products of single-qubit 
predictions. 

The inconsistency between joint and correlated local 
predictions is a general feature of our LHV model for n- 
qubit GHZ states. It occurs only for joint measurements 
that involve Y measurements on some qubits and that 
have a definite outcome, i.e., joint measurements of sta- 
bilizer elements that contain y's. Joint measurements 
that yield a random result do not suffer from this prob- 
lem because the randomness of a product is unaffected 
by discarding i's. More precisely then, the inconsistency 



occurs only for joint measurements that are products of 
X's and y's on all the qubits, with the number y's being 
an even number that is not a múltiple of 4. 

The protocol for ensuring consistency between joint 
and composite local predictions proceeds as follows. An 
observer called Alice, stationed at, say, the first qubit, is 
put in charge of ensuring consistency with single-qubit 
measurements. She does so by changing or not chang- 
ing the sign of the outcome on her qubit, based on what 
is measured on her qubit and information she receives 
about what is measured on the other qubits. Because 
of the qubit-exchange symmetry of the GHZ state (and 
of the LHV table), an observer stationed at any qubit 
could play the role of Alice. Alice ensures consistency by 
changing the sign of her local outcome if and only if (i) a 
measurement of X or Y is made on her qubit and (ii) the 
total number of Y measurements on all qubits is an even 
number that is not a múltiple of 4. The protocol requires 
n — 1 bits of communication as each of the other qubits 
reports to Alice whether Y was measured on that qubit. 
The protocol clearly fixes all those cases that need cor- 
rection; just as important, in all situations where Alice 
flips her qubit, all subsets of qubits that include Alice's 
qubit, except for the case of a needed correction, have a 
random measurement product, which is therefore unaf- 
fected by Alice's flip. The success of this protocol clearly 
relies on very special properties of the stabilizer group 
for the n-qubit GHZ state. 

We can put the protocol in a more mathematical form 
by letting r± = 1 if an X or Y measurement is made on 
the first qubit and r\ — otherwise and by letting qj = i 
if y is measured on the jth qubit and qj = 1 otherwise. 
Alice ensures consistency by flipping her local outcome if 
and only if p n — r\q± ■ ■ ■ q n = — 1. This formulation al- 
lows us to see easily that we can do a bit better than the 
n — 1 bits required by the original protocol. The key is to 
notice that when p n — ±i, all subsets of qubits that in- 
clude Alice's qubit have a random measurement product, 
so a flip by Alice goes unnoticed. As a result, Alice can 
get by with the truncated product p„_i = r\q\ ■ ■ ■ q n -\, 
flipping her local outcome if and only if p n -\ = i or — 1. 
This scheme requires the promised n — 2 bits of commu- 
nication, because Alice doesn't need to know whether a 
y measurement is made on the nth qubit; it works be- 
cause Alice flips whenever p n = —1, as required, with the 
additional flips when p n = ±i not doing any harm. 

The consistency scheme generalizes trivially to the case 
of Pauli-product measurements made on l disjoint sets 
of qubits. For each set k chosen from the l sets, the 
table yields a measurement product that is the predicted 
outcome multiplied by qu = i or q% = 1. Putting Alice 
in charge of the first set, all but the last of the other sets 
communicates qk to Alice, who computes the product 
?"i<Zi • • " <7í-i) where r\ = if no measurement or a Z 
measurement is made on any qubit in her set and r\ = 1 
otherwise. Alice flips her set's outcome if and only if 
T\q\ • • • qi—i = i or — 1. Consistency is thus ensured at 
the price of l — 2 bits of communication. 
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IV. PROOF OF ÒPTIM ALITY 

Using an elaboration of Mermin's GHZ argument , 
we now demonstrate that our model is optimal by show- 
ing that any protocol that is allowed at most n — 3 
bits of classical communication is incapable of yielding 
all quantum-mechanical predictions for measurements of 
Pauli products on \ip n ) and their correlations. For this 
purpose, imagine the n qubits as the nodes of a graph; 
two qubits are connected by a line if at least one bit is 
communicated between them. The graph partitions the 
qubits into disjoint connected subsets. There being at 
most n — 3 lines, it follows that there are at least three 
disconnected subsets, since at best each line consolidates 
two subsets into one, thcreby eliminating one subset. 
Moreover, it is always possible to arrange the commu- 
nication so that there are three subsets. The communi- 
cation can do no more than allow us to treat all Pauli 
products within a subset as a single joint observable. We 
can restrict attention to the case of three subsets, since 
amalgamating disconnected subsets allows the communi- 
cation more power than it actually has. 

The situation then is that we have three disconnected 
subsets, containing k, l, and m qubits, with k + l+rn = n. 
Since the GHZ state is invariant under qubit exchange, 
we can make the first k qubits those in the first sub- 
set and the next l qubits those in the second subset, 
leaving the final m qubits to be those in the third sub- 
set. We now define six Pauli products: A = X^ k ^X 
and B = X^ k -^Y for the first subset, C = X^-^X 
and D = X^^^^Y for the second subset, and E = 
X®(m-i)x and F = X^ (m -^Y for the third subset. 
The four operators, ACE, —ADF, —BCF, —BDE, are 
in the stabilizer group of \ip n ) and thus give a definite 
outcome +1 when measured, implying that a measure- 
ment of any two of the operators in the product can be 
used to predict with certainty the result of a measure of 
the third. The EPR reality criterion then says that we 
should associate elements of reality, having vàlues ±1, 
with the six Pauli products A-F. Denoting the vàlues 
of these elements of reality in the obvious way, the def- 
inite vàlues of the last three stabilizer elements imply 
that adf = bcf = bde = — 1. The product of these three 
quantities is ace = — 1, contradicting the +1 prediction 
of quantum mechanics for a measurement of ACE. 

For complcteness, we note another form of the argu- 
ment. According to the elements of reality, the observ- 
able M = ACE - ADF - BCF - BDE has the value 

Ai = ace — adf — bcf — bde = c(ae — bf) — d(af + be) . (1) 

Since ae — ±bf af = ±òe, it is easy to see that 
Ai = ±2. This implies that the expectation value sat- 
isfies |(.M)| < 2, whereas the n-qubit GHZ state has 
(Ai) = 4. This form of the argument does not make 
use of the properties of the GHZ state, and it makes 
clear that stochastic models can do no better than the 
deterministic models considered here. We also note that 
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FIG. 1: Quantum circuit that generates the three-qubit GHZ 
state. 



this argument pro duces a Bell inequality with auxiliary 
communication 16] . 



V. QUANTUM CIRCUIT 

Tablc ^ is the basis of our communication-assisted 
LHV model. It arises naturally from a quantum circuit 
that creates the n-qubit GHZ state from an initial state 
1 00 . . . 0) . One such circuit consists of a Hadamard gate 
on the first qubit followed by n — 1 C-NOT gates, with 
the leading qubit being the control and the remaining 
qubits serving successively as targets. The three-qubit 
version of this circuit is shown in Fig. The Hadamard 
gate H transforms the Pauli operators according to 



HXHÏ = Z, HYW = -Y, HZH 1 = X 



t _ 



Similarly, under the action of C-NOT, we have 



(2) 



C(XI)CÏ = XX , C(FJ)Ct = YX , C{ZF)CÏ = ZI , 

c(ix)c f = ix , c*(/y)c t = zy , c(iz)c f = zz , 

(3) 

where the first qubit is the control and the second is the 
target. These operator transformations lead to the table 
update rules given in Fig. [3 which traces the evolution of 
the LHV table during the creation of a three-qubit GHZ 
state using the circuit of Fig. ^ A simple generalization 
to n qubits leads to Table ^ for the n-qubit GHZ state. 

The C-NOT update rules given in Fig. [21 must be 
consistent with the fifteen nontrivial transformations of 
Pauli products generated by C. Six of these transfor- 
mations, listed in Eq. ©, serve as the basis for the up- 
date rules. Because C — C* , the rules are automatically 
consistent with four other transformations. In addition, 
the rules are clearly consistent with the transformation 
C(ZX)C^ = ZX. Consistency with the remaining four 
transformations, C (XY) = YZ, C (XZ) C^ = -YY , 
and their inverses, requires that 



v^v^ vF yF -w-I yl yl 

■ A C ï t — r c — ï C Z/ C Zj t jL t ! 

X\Z\ = -YjY? = -YlZ\X\Y} 



(4) 



These relations do not hold generally, but are satisficd 
if the initial entries for both the control and target are 
correlated according to XY Z = i (or XY Z — —i), with 
X and Z real and Y imaginary. These conditions hold in 
all our applications of C-NOT. It is for this reason that 
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X Y Z X Y Z X Y Z 

i?i — iRi 1 Hi 1 iRi R 1 C-NOT12 R2 1R1R2 Ri 

R 2 iR 2 1 * Ri ÍR2 1 ' R2 ÍR1R2 Ri 

R 3 iR 3 1 R 3 iR 3 1 R 3 iR 3 1 

X Y Z 

C-NOT13 -R2-R3 Í-R1-R2-R3 Ri 

R2 ÍR1R2 Ri 

R 3 iRiR 3 Ri 

FIG. 2: Evolution of the LHV table during the creation 
of a three-qubit GHZ state using the circuit of Fig. Q 
The initial table yields the correct quantum predictions for 
the state |000). The rules for updating the table through 
Hadamard and C-NOT gates come from the operator trans- 
formations @ and J3J|. The Hadamard update rules are 

x F = z 1 , y f = -y 1 , z F = x 1 , 

where I and F denote the initial and final vàlues of a table en- 
try, before and after the application of the gate. The rules for 
updating through a C-NOT, with control c and target í, are 

Y F — Y 1 Y 1 V F — V 1 Y 1 7 F — 7 1 
yF _ y 1 V F — 7W 1 7 F — 7 l 7 l 

The update rules are local in that they only require changes to 
table entries corresponding to the qubits involved in a gate. 
The subscripts on the random variables in Table Q are now 
seen to represent the qubits to which these variables were ini- 
tially associated. Correlations arising from the application of 
C-NOT gates then correspond to pairs of identical subscripts. 

the initial sign of the Y entry for the fixst qubit (see the 
first table in Fig. 0) is opposite that of the remaining 
qubits. 

VI. CONCLUSION 

We have shown that it is possible to reproduce cor- 
rectly the quantum-mechanical measurement predictions 
for the set of all n-fold products of Pauli operators on an 
71-qubit GHZ state using only Mermin-type LHV's and 
n — 2 bits of classical communication. The n — 2 bits of 
communication, shown here to be optimal, are required 
to ensure that the products of local measurement pre- 
dictions are consistent with the corresponding joint pre- 
dictions. We also show how our model arises naturally 



from a simulation of a quantum circuit that creates the 
?i-qubit GHZ state. 

The circuit that creates the GHZ state is an examplc 
of the Gottesman-Knill circuits mentioned in the Intro- 
duction. Although GK circuits can produce global en- 
tanglement, as in the n-qubit GHZ state, they can be 
emciently simulated in 0(n 2 /logn) steps on a classical 
computer [3, 1 1 5| . The important difference between 
the GK simulation of the circuit that creates the n-qubit 
GHZ state and our communication-assisted LHV simu- 
lation is that the GK algorithm tracks the evolution of 
nonlocal hidden variables, represented by the n gener- 
ators of the stabilizer group, and therefore requires no 
communication overhead to give correct predictions. 

A generalization of our results might lead to a new per- 
spective on the GK theorem since our results imply that, 
at least in this limited case, we can replace the nonlocal 
hidden variables represented by the stabilizer generators 
with LHV's and a linear amount of classical communi- 
cation. We conjecture that this is a genèric feature of 
quantum circuits obeying the constraints of the GK theo- 
rem; that is, we expect that any quantum state produccd 
by a GK circuit can be modelled with EPR elements of 
reality plus an amount of classical communication that 
scales linearly in the number of qubits. 

Thcre are two main obstacles to a straightforward ex- 
tension of our LHV model to general GK circuits. The 
first is the difficulty, expressed in Eq. Q , in maintaining 
the consistency conditions for the C-NOT update rules. 
The second is the reliance of our communication protocol 
on very special properties of the GHZ state; for general 
GK states, the communication protocol will have to be 
more complicated, with a proof of optimality correspond- 
ingly more complicated as well. Nonctheless, the exis- 
tence of a communication-assisted LHV model for arbi- 
trary GK circuits and the entangled states they produce 
is currently under investigation. 
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